Andreev quantum dots for spin manipulation 
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We investigate the feasibility of manipulating individual spin in a superconducting junction where 
Bogolyubov quasiparticles can be trapped in discrete Andreev levels. We call this system Andreev 
Quantum Dot (AQD) to be contrasted with a common semiconductor quantum dot. We show that 
AQD can be brought into a spin- 1/2 state. The coupling between spin and superconducting current 
facilitate manipulation and measurement of this state. We demonstrate that one can operate two 
inductively coupled AQD's as a XOR gate, this enables quantum computing applications. 

PACS numbers: 05.60.Gg,03.67.-a,73.21.-b 
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Manipulation and operation of individual quantum 
systems and arrays of such systems, so-called "quantum 
machines" is now in focus of both experimental and theo- 
retical research Jl Q . The progress in quantum comput- 
ing algorithms |3j has demonstrated potential applicabil- 
ity of quantum mechanics thus stimulating various pro- 
posals to implement arrays of operational two-state sys- 
tems (qubits) in solid-state 0, 1111,1111 0E1 El El Gj 
Many proposals concern quantum dots. The quantum 
dots are often referred to as artificial atoms since they 
confine a discrete number of particles that occupy dis- 
crete quantum states. In contrast to atoms, the proper- 
ties of quantum dots can be tuned and their charge and 
spin degrees of freedom can be controlled. This would al- 
low for quantum manipulation. An interesting and elabo- 
rated proposal 4, ;5| utilizes spin states of semiconductor 
quantum dots. However, the complexity of the manipu- 
lation schemes proposed and severe difficulties with the 
read-out of these spin states [1 Q drives one to think of 
alternatives. 

Below we present an alternative scheme for individ- 
ual spin manipulation. We concentrate on sufficiently 
resistive superconducting constrictions where individual 
Bogolybov quasiparticles can be trapped in discrete An- 
dreev bound states. We refer to such system as Andreev 
Quantum Dots (AQD). An AQD resembles a common 
quantum dot as long as discreetness of a (quaisi)particle 
number, spectrum and spin is concerned. Albeit in con- 
trast to a common quantum dot the charge of the AQD 
is not fixed. This allows for superconducting current in 
the constriction and makes electron-electron interaction 
negligible. 

We propose to utilize spin states of the AQD's. We 
show that an AQD can be brought to the state with spin- 
1/2 that persist over a long time. It is important that the 
spin direction in this state determines the superconduct- 
ing current in the constriction, thus solving the read-out 
problem. We demonstrate that the spin state of a sin- 
gle AQD can be manipulated. Further, the two dots can 
be inductively coupled to make a XOR quantum gate. 



Quantum information theory Q proves that this enables 
one to build a universal quantum computer. 

The AQD can be formed in any constriction between 
two superconducting leads that have a gap in energy 
spectrum. If an electron with the energy below the 
gap tries to escape to the bulk of a superconductor, 
it is reflected back as a hole (Andreev reflection [lfij|). 
which also can not escape due to the same reason. So 
that, the junction confines quasiparticles that are co- 
herent mixtures of electron and holes. Their discrete 
energy levels and eigenfunctions are determined by Bo- 
goliubov equations (BdG) Ell- It i s sometimes forgot- 
ten that these equations do possess a spin structure. 
Bogoliubov eigenfunctions are made of two spinors 17| 
u a , v a [coefficients of the Bogoliubov transformation lfj 
*(r,o-) = E„K(r,cr) 7 „ + gv^v^r, that satisfy 



eu a = H a u 



At 



ev a = -[H*] 



(1) 
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Here "hat" denotes an operator over orbital degrees 
of freedom. We make explicit the spin structure of 
the single-particle Hamiltonian H and pair potential A, 
g a p = ia y being metric tensor in spinor space E3, 
{H*) a p = g va {H v ^g^. By virtue of Eq. quasipar- 
ticle energy levels always come in pairs: each eigenstate 
with energy e has a counterpart with energy — e. This 
is due to a double-counting: there are two quasiparticle 
eigenfunctions per each state of H . Should H possess no 
spin structure, Andreev levels are spin-degenerate. For 
many problems that do not involve spin, one can avoid 
the double counting by considering one non- degenerate 
level. This technical trick does not corres pon d to original 
formulation of superconductivity theory |lo . Il8j , neither 
it gives the correct description of spin in superconductors. 

We concentrate on a short constriction, such that the 
typical time for an electron r (li M to traverse the junc- 
tion satisfies the condition T flight <C h/A. In the limit 
r m , t A/H — > 0, and in the absence of magnetic field An- 
dreev levels are spin-degenerate and can be universally 
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expressed through eigenvalues T n of the transmission 

matrix square, e„ ljn2;CT = sign(n 2 )Ay / 1 - T ni sin 2 (<^/2). 
Here the integer index n\ labels orbital channels, n-i — 
±1, a = ±1 is spin- index and ip stands for the supercon- 
ducting phase difference between the leads. 

Andreev levels that are relevant for electron trans- 
port, and for manipulation of spin states, originate from 
T n < 1. These levels are distributed in energy strip 
A| cos(iy9/2)| < |e| < A. Their typical spacing is given 
by SE ~ AGqR, R being the normal state resistance of 
the constriction, Gq being the conductance quantum [2(| • 
This estimation is valid for any sufficiently disordered 
constriction except tunnel junctions for which all T n <C 1. 
To give an estimation, for R f« 10 Ohm and A « 10K, 
the spacing SE w lOOmK, this is typical for semiconduc- 
tor QD. ft is remarkable that that the Andreev levels do 
not depend on microscopic and geometric details of the 
constriction by means other than conductance and trans- 
mission amplitudes. This fact considerably simplifies the 
fabrication of Andreev dots with discrete Andreev levels. 

In the ground state of the dot, quasiparticles oc- 
cupy Andreev levels with negative energy. The in- 
dependent part of the ground state energy reads E a — 
1/2 J2na £ no-Q(— £na) where we sum over channel and 
spin index; O(x) = 1 for x > and zero for x < 
0. The factor 1/2 comes from the double-counting 
mentioned. 01 If spin splitting of Andreev levels is 
smaller than SE, the ground state has zero spin, since 
both components of the spin dublet are occupied by 
quasiparticles. In an excited state of the AQD, some 
Andreev levels with positive energy are populated. Let 
us concentrate on a given transport channel where there 
are two such Andreev levels corresponding to two spin di- 
rections. One quasiparticle fills either level, the AQD has 
spin- 1/2. Second quasiparticle fills the level with the op- 
posite spin resulting in an excited spin-singlet state. (Fig. 
la) The spin-1/2 state of an AQD with the lowest en- 
ergy (that corresponds to the most transparent transport 
channel) is of particular interest because it is very stable. 
The transition to ground state require the 1/2 change of 
spin. This means that a quasiparticle must either leave 
or enter the AQD. The probabilities of these processes 
contain exponentially small factors exp(— A/fceT), this 
means that at zero temperature the AQD would remain 
in spin-1/2 forever. The physics involved is very simi- 
lar to well-known parity effect in superconducting grains 
[2l| . Thus it is possible to preserve the system for a long 
time in this spin-1/2 state. 

How to set the AQD to spin-1/2 state? Possibili- 
ties include microwave irradiation [23] and quasiparticle 
injection [23)]. We concentrate on the first possibility. 
Let us assume that the irradiation frequency to satisfies 
A + Eq < fttjj < 2A. Under these conditions, the ab- 
sorption of the irradiation quanta takes place in the con- 
striction only. If initial state is the ground state, each 
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FIG. 1: a) Excited states of Andreev quantum dot in a given 
transport channel, b) The energies of AQD states correspond- 
ing to a given transport level versus phase difference. The 
lower and upper curve correspond to ground state and ex- 
cited singlet respectively. The middle curves correspond to 
two spin-1/2 excited states, spin splitting being due to spin- 
orbit interaction. 



absorption generates two quasiparticles. There are two 
possibilities: 1) both quasiparticles appear in bound An- 
dreev states; 2) one quasiparticle appears in a bound 
state whereas another one acquires energy > A and gets 
to the extended state (photoemission); the latter quasi- 
particle leaves the AQD almost immediately and never 
comes back. Next absorption process can occur in ex- 
cited state, in this case the energy balance allows extra 
processes: 1) two quasiparticles leave the dot, 2) one 
quasiparticle in a bound state is excited over the energy 
barrier and leaves the dot. In any case, the processes of 
the first type are irrelevant not changing the parity of 
the AQD. Processes of the second type switch the AQD 
between states with even and odd number of quasiparti- 
cles. If the irradiation lasts long enough for many pro- 
cesses of the second type to occur, the AQD is in the 
state with half-integer spin with 50 % probability. Let 
us now switch off the irradiation. If there is an even 
number of quasiparticles in the AQD, the subsequent en- 
ergy relaxation will drive the system to the ground state. 
For an odd number of particles, the relaxation will result 
in a stable single quasiparticle occupying the lowermost 
Andreev level, the lowermost spin-1/2 state we are after. 

How to detect the spin-1/2 state? It is important 
that the superconducting current in the constriction is 
different for the states involved, since their ^-dependent 
energies arc different. Here we concentare on a chan- 
nel n and count energy from the ground state energy of 



quasiparticles in all other channels, E g 



The en- 



ergies then are — |e n |, 0, and |e„| for ground, spin-1/2, 
and excited singlet state respectively. (Fig. lb) and the 
superconducting current equals to J = ed v E/h. The 
change from the ground to spin-1/2 state is therefore 
manifested as a change of supeconducting current by a 
value of 51 = ed^e na /H. The detection of such current 
jumps in superconducting constrictions would amount to 
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the direct experimental observation of the spin-1/2 state. 

How to detect spin in the spin-1/2 state? An impor- 
tant advantage of an AQD is that its spin state affects 
the superconducting current, the latter being detected 
and/or measured. This is due to spin-orbit splitting of 
Andreev levels. Generally, one expects such splitting 
given the fact that the supeconducting phase difference 
changes sign under time reversal so that Kramers the- 
orem does not hold. A confusing circumstance is that 
the Kramers theorem does hold in the universal limit of 
short constriction considered. To get spin-orbit splitting, 
one considers extra perturbative corrections of the first 
order in r~. ..A/H. The calculation based on the scat- 

flight / 

tering matrix approach 19] with the use of Eq.QJ yields 
the following effective Hamiltonian for the spin-spitting 
in question: 

E n SO) = A ( a n ■ tr) sin(y)(r flight A/?j), (2) 

a being the pseudovector operator of spin. Here a di- 
mensionless pseudovector a n is a property of a given 
Andreev level not depending on tp; it is proportional 
to the spin-orbit constant that routinely contains a nu- 
clear charge of the material Z, \a\ n Z(e 2 /hc). For ex- 
ample, in a quasi-ballistic SNS junction a is directed 
along the vector product of the quasiparticle momentum 
P (IpI ~ Pf) and the direction of the electric (crystal) 
field E for tp > 0: a J] p x E; but when tp < 0, p 
changes its direction to the opposite and a |f p x E. 
This change of the direction of the quasiparticle momen- 
tum (Andreev state "chirality" ) with the sign of ip is the 
reason why the Andreev level spin-splitting in Eq. J2J) has 
rather unusual odd dependence on tp [this is general prop- 
erty of spin-orbit Andreeev level splitting]. The techni- 
cal model to derive Eq.(|2J) was a one-channel conductor 
with two scatterers separated by a distance correspond- 
ing to Tflight- Both scattering matrices contained spin- 
orbit part. The model showed that in not very short 
constrictions (T !Ught A/h ~ 1) spin-orbit level splitting 
aquired more complicated dependence on tp and could 
be of the order of A provided spin-orbit scattering was 
comparable with orbital one. Thus the splitting can be 
huge: for a material with heavy nuclei a can become of 
the order of unity. To show how big can be the effect, let 
us consider an AQD embedded into a superconducting 
loop with self- inductance L. The spin produces an extra 
superconducting current I a — ead v E^ so ^ h and an ex- 
tra magnetic moment per spin thereby. Let us estimate 
the maximum possible value of this magnetic moment. 
For this, we set r fligM ~ A/h ,\a\ ~ 1, and concentrate 
on a resistive constriction R i2 Rq [so that the critical 
current is of the order of e A / h] . The inductance L of a 
typical SQUID loop does not exceed 1cm [3] thus the 
extra magnetic moment SM < (eA/K)L 2 / c ~ 10 12 /j,b- 

Let us give a simple example of spin manipulation 
in the AQD. For this, we need to invoke Zeeman split- 
ting of Andreev level in magnetic field aE^ z > in addition 



to spin-orbit splitting aE n ■ This allows us to control 
the direction of spin quantization by the magnetic field, 
this possibility being absent if we work with spin-orbit 
only. To achieve comparable E^ z ' and E^- so \ the mag- 
netic field should be almost "in plane" not affecting the 
flux in the SQUID loop and thus tp. Let us assume that 
we manage to achieve this so that the quantization axis 
of spin may deviate substantially from a. The spin wave 
function is therefore a coherent mixture of the states | f), 
| J.) with spin parallel or antiparallel to a respectively. 
Quite generally, * = a\ f) + b\ j), where |a| 2 + |6| 2 = 1. 
If there is no in-plane field, the system is in the ground 
state | J,), spin quantization axis being parallel to a. Now 
let us switch on the in-plane field. The Hamiltonian gov- 
erning dynamics of the wave function will thus become : 
H = guB((T-H)+E(go)- Let us assume the simplest form 
of the resulting Hamiltonian: H = E^a y + E^ so ^a z , 
E ( z ) ~ E( so l The wave function will then evolve ac- 
cording to 

= | J.) (cos(fif) + isin(fit)(^ so Vfifi)) - 

-| ■\)sin{flt){E {z) /hn), 

where hfl = ^/(E^ so )) 2 + (E^) 2 is the frequency of 
Rabi oscillations. It is important to note that these os- 
cillations can be readily detected since they produce an 
alternating current 

I a (t) = 2 sin 2 (M) (^pj ed v EV°)/h. (3) 

Another way of manipulation is readily borrowed from 
the quantum optics: if the in-plane magnetic field oscil- 
lates with the resonant frequency Ego , & significant ma- 
nipulaton effect can be achieved even if E^ < 
[25j. A general unitary transformation of the spin- wave 
function can be performed exposing the junction to the 
time-dependent E^ z ', quite similar to many other solid- 
state implementations of the qubits. 

Let us discuss now how the Andreev quantum dots can 
be utilized for universal quantum computations. An 
AQD in the spin-1/2 state would be a qubit. We have dis- 
cussed above how to manipulate the spin of a single AQD. 
This is how the single-qubit operations can be performed. 
A quantum computation algorithm should involve two- 
qubit operations as well. An important theoretical result 
3] establishes that the XOR operation with two qubits 
along with single-qubit operations forms a minimum set 
that is sufficient to build up an arbitrary complex uni- 
versal quantum computer. The XOR operation does the 
following: given two qubits in the states |a;), \y), it leaves 
the \y) state unchanged if |x) = | f), while flipping it 
when \x) = \ |). So we concentrate on a possible realiza- 
tion of XOR operation for two AQD's. 

The basic idea is to organize the interaction between 
AQD's via inductive coupling between SQUID loops con- 
taining these AQD's [see Fig. [2]. The two-qubit opera- 
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FIG. 2: The interaction between AQD qubits is organized by 
means of variable inductive coupling. This allows for further 
arraying. The circuit shown can be used as an XOR gate. 



tions are performed by varying the mutual inductance 
between two given loops. Indeed, in this case the inter- 
action between two AQD's can be described by a simple 
hamiltonian H — Lnlxli, £12 being the mutual induc- 
tance, ^i(2) standing for operators of currents in corre- 
sponding loops. Owing to spin-orbit interaction, each 
current depends on spin state of corresponding AQD: 

Ii t 2 = /i 2+/{ 2 s i 2- Here we choose z-axis in the space of 
each spin in the direction of corresponding pseudovector 
a. This is to stress the following circumstance: although 
the pseudovectors a in different AQD's may differ, there 
is always only one component of spin that is reflected in 
the current and therefore takes part in the interaction. So 
that, the relevant part of the Hamiltonian can be written 
as 



H = H\s\ + H2S2 + H\is\s\. 



(4) 



This simple Ising-type form of Hamiltonian brings us 
to the old-fashioned but solid "optical" quantum com- 
puter |26j |. In this approach, the one-bit operations are 
performed at H\i = by pulses at resonant frequen- 
cies Hi/h or H2/h, the pulse duration being tuned to 
achieve the spin flip. The XOR operation is performed at 
H12 ^ by the same pulse with frequency (Hi + B.2)l%. 
An alternative way is to use non-oscillating pulses of 
H12 1 h. Such pulses would shift phases of two states with 
antiparallel spins with respect to the phases of the states 
with parallel spins thus realizing "quantum phase shift 
gate" |13 ■ The XOR operation can be performed by 
combining two such phase shifts with two rotations of 
the target spin. 

This approach of organizing two-qubit interactions has 
two important practical advantages. First, in contrast to 
other spin-based proposals, the interaction does not have 
to be organized at microscopic level. To exaggerate, one 
can use inch-scale transformers to vary inductive coupling 
between the AQD's. To make a practical suggestion, one 
can use the well-developed techniques of SQUID circuitry 
[28L Eif to couple, array, bias, and measure many AQD 
qubits. Second advantage is the simple Ising form of the 



resulting interaction that prevents undesired phase shifts 
and simplifies design of complicated quantum circuits. 

To conclude, we analyze prospectives of Andreev quan- 
tum dots for spin manipulation and quantum computing. 
Our theoretical results seem to be promising enough to 
launch detailed experimental investigations and design 
efforts in this direction. 
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